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ABSTRACT 

We show that there are strong local correlations between metallicity, surface brightness, and dynam- 
ical mas s-to-light ratio within M33, analogous to the fundamental line of dwarf galaxies identified by 
iPrada fc Burkert ( 2002). Using near-infrared imaging from 2MASS, the published rotation curve of 
M33, and literature measurements of the metallicities of H II regions and supergiant stars, we demon- 
strate that these correlations hold for points at radial distances between 140 pc and 6.2 kpc from the 
center of the galaxy. At a given metallicity or surface brightness, M33 has a mass-to-light ratio ap- 
proximately four times as large as the Local Group dwarf galaxies; other than this constant offset, we 
see broad agreement between the M33 and dwarf galaxy data. We use analytical arguments to show 
that at least two of the three fundamental line correlations are basic properties of disk galaxies that 
can be derived from very general assumptions. We investigate the effect of supernova feedback on the 
fundamental line with numerical models and conclude that while feedback clearly controls the scatter 
in the fundamental line, it is not needed to create the fundamental line itself, in agreement with our 
analytical calculations. We also compare the M33 data with measurements of a simulated disk galaxy, 
finding that the simulation reproduces the trends in the data correctly and matches the fundamental 
line, although the metallicity of the simulated galaxy is too high, and the surface brightness is lower 
than that of M33. 

Subject headings: galaxies: abundances — galaxies: formation — galaxies: fundamental parameters 
— galaxies: individual (M33) — galaxies: kinematics and dynamics — galaxies: 
spiral 



1. INTRODUCTION 

'Pr ada fc Burkerti lj2002l hereafter PB02) recently 
found that the dwarf galaxies in the Local Group 
obey a puzzling correlation between mean metallicity 
and global dynamical mass-to- light ratio (M/L). By 
combining this relationship with the previously known 
tren d of metallicity with central surface brightness 
(e.g, 'Phi UipDS. Edmunds, fc Da^lToM iCaldwell et all 
^98), PB02 demonstrated that all three parameters to- 
gether form an even stronger relationship, which they 
labeled the fundamental line of dwarf galaxies. PB02 
showed that the metallicity-M/L ratio correlation can be 
explained by a simple chemical enrichment model. In this 
model, galactic winds continuously remove the heavy el- 
ements formed in dwarf galaxies from their interstellar 
medium, and star formation in these galaxies ends when 
all of the gas has been blown out of the system. Be- 
cause more metals are retained by more massive galax- 
ies (which have lower mass-to-light ratios) this model is 
able to reproduce the observed metallicity-M/L relation- 
ship. However, the model does not make clear why in- 
cluding surface brightness as a third parameter improves 
the correlation, and it is also surprising that variables 
such as morphological type and star formation history 
have no apparent effect on the fundamental line. PBG2 
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proposed that feedback from star formation or super- 
novae could be involved in regulating the fundamental 
hue, and Dekcl & Woo (2003) and Tassis ct al. (2003) 
used analytical arguments and simulations, respectively, 
to show that supernova feedback can produce relation- 
ships like the fundamental line, but additional work is 
needed to improve this understanding. 

Two natural questions arise from these previous inves- 
tigations of the fundamental line. First, does the funda- 
mental line hold for massive galaxies as well as dwarfs? 
Galactic winds and supernova feedback should operate 
less effectively in larger galaxies as a result of their deeper 
potential wells, so if either of these processes is the pri- 
mary mechanism responsible for setting the fundamen- 
tal line correlations, the fundamental line should not be 
the same for massive systems. Second, is the fundamen- 
tal line a local law as well as a global one? If not, it 
would be puzzling that all dwarf galaxies know to fol- 
low the relationship globally without having any corre- 
sponding correlations on smaller scales. The data to test 
the fundamental line locally in dwarfs do not yet ex- 
ist, because spatially resolved metallicity measurements 
have not been obtained, and few dSph and dE galaxies 
have well-resolved kinematics either. Answering either of 
these questions with existing data therefore requires us 
to move to more massive galaxies. 

In order to explore further the nature of the funda- 
mental line, we consider the case of a well-resolved disk 
galaxy with known variations of its metallicity, surface 
brightness, and mass-to-light ratio with radius. The best 
example of a galaxy with these necessary characteristics 
is M33. We will address the following questions with 
this study. 1) Does the fundamental line hold locally, at 
every point within a galaxy? 2) Do the parameters of 
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the global and local (assuming one exists) fundamental 
lines agree? 3) What is the origin of the fundamental 
line? 4) Can simulations of galaxy formation reproduce 
the fundamental line? 

In this paper, we use several recently published data 
sets to investigate the local fundamental line in M33. In 
the next section, we combine a dee p Two Micr on All 
Sky Survey (2MASS) mosaic of M33 (IBlock et al.ll200l 
high resolution CO and H I rotation curves (ICorbellil 
l2003HCorbelh fc Sa luccr2000l and literature metaUicity 
measurements to construct the necessary observational 
quantities. In Sj^we perform fits to the various relation- 
ships between mass-to-light ratio, metallicity, and sur- 
face brightness, and in 21 '^fi compare the locally deter- 
mined fundamental line of disk galaxies to the fundamen- 
tal line of dwarf galaxies, give derivations of two of the 
fundamental line relationships, and discuss the implica- 
tions of our findings. We also use the r ecent ACDM disk 
galaxy simulation bv lRobertson et aP l)2004D to test how 
accurately current simulations can reproduce the prop- 
erties of galaxies and we employ additional toy models 
of galaxy formation to understand the impact of feed- 
back on the fundamental line. We briefly summarize our 
results and conclude in ^ 

2. DATA AND ANALYSIS 

Of the observations required for this study, the most 
restrictive is the large sample of accurate metallicity mea- 
surements as a function of radius. Many nearby galax- 
ies have high-quality rotation curves and deep optical 
imaging, although Ks band imaging to the depth of the 
2MASS images (see W2.2^ is less common. However, only 
two galaxies — MlOl IjKennicutt. Bresolin. fc GarnettI 
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|2003|) and M33 — have measured abundances for more 
than ~ 10 H II regions. Because MlOl is a massive 
galaxy (and therefore more dynamically dominated by 
baryons), is nearly face-on (making its kinematics more 
difficult to measure accurately), and is significantly lop- 
sided (indicating a possible recent interaction), it is not 
an ideal target for this study, and we instead focus ex- 
clusively on M33. 

2.1. Abundances 

Studying the fundamental line over a large range of 
galactic environments requires abundance measurements 
covering as wide a range of radii as possible. Because we 
will be comparing abundances to other local properties, 
it is important that we use a consistent abundance scale. 
We therefore limit our H II region sample to objects with 
directly-determined electron temperatures (although we 
do include CC93, the closest H II region to the center 
of M33, for which abun dances were deri ved from strong- 
line model fitting bv iVflchez et al.lll988(l . We also utilize 
abundances for M33 B supergiants, which have metal- 
licities consistent with those of the H II regions in M33, 
although with slightly larger scatter. 

The abundance of oxygen in the interstellar medium 
(ISM) and in the young stellar population of M33 
has been sampled across the disk of the galaxy using 
both star-forming regions and young stars. We have 
compiled abundance data for 13 H II regions derived 
from optical spectroscopy l) Smith ' 1975; Kwitter & Alleil 
Il98lt IVflchez et aLllT988l) and ISO LWS spectroscopy 
(jHiedon et alJl2003j) . The compilation of abundances for 
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Fig. 1. — Oxygen abundances as a function of radius in M33. The 
filled points represent abundance measurements from H II regions 
and the open points represent abundance measurements from B 
supergiants. 



13 bri g ht M33 sup e rgiant s cam e from iMonteverde et al] 
l|T997|) . lUrbaneia l l|200^ . and lUrbaneia et al. I l|200fi() . 

These works used quantitative spectroscopy of M33 su- 
pergiants, subsequently analysed with state of the art 
model atmospheres of massive stars. Combining these 
data sets yields 26 metallicity measurements in M33, 
spanning radii from 140 pc out to 6.2 kpc. The metallic- 
ity measurements employed in our study are plotted as a 
function of radius in Figure^ Note that we assume a dis- 
tance of 800 kpc througho ut this paper (Lee__et al. 20Q2 
iMcConnachie et alJl200i|) . although some other recent 
determinations suggest a distance of in ore than 900 kpc 
l|Kim et alJl200l ICiaTdullo et alJl200^ . 

Abundance data for other astrophysical sources have 
been reported in the literature for M33: e^., plane- 
tary nebulae (PN; Magrini et al. 2003, 20o4_and su- 
pernova remnants fSNR,: Blair fc Kirshnenil985^ : overall 
the PN and SNR abundances show a general agreement 
with the abundances derived for the H II regions (see 
[Vflchez et al. 1988; Magrini ct al. 2004; Stasiilska et aO 
|2005|) . However, only small samples of these objects have 
been observed, and the scatter in the data is large, in- 
dicating that further observations are needed to under- 
stand how accurately they trace the present ISM abun- 
dances. Therefore, for the purposes of this paper we have 
not used PN and SNR data to trace the abundance gra- 
dient, though we point out that PN and (possibly also) 
SNR abundance results appear consistent with the data 
set used here. 

2.2. 2MASS Images 

In addition to providing relatively shallow near- 
infrared images of the entire sky, the 2MASS project 
also made deeper observations of several nearby galax- 
ies. The deep exposures of M33 have integration times 
six times longer than the standard 2MASS data, so they 
reach ~ 1 mag fainter in surface brightness (" Block et alJ 
120041) . Tom Jarrett has kindly provided us with a mosaic 
of the individual images that covers the entire galaxy, 
with foreground stars and background (sky) variations 
removed. 
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Ks band surface brightnesses at the position of each 
H II region and B star were measured by adding up the 
flux contained in circular apertures of radius 7". The 
local surface brightnesses measured in this way are con- 
sistent with the azimuthally averaged surface brightness 
(see below) at the same radius within the uncertainties. 
Although the H II regions themselves contribute negligi- 
bly to the observed flux at this wavelength, in some cases 
a single star coincidentally located very close to the H II 
regions adds significantly to the flux, biasing the surface 
brightness high. To avoid this contamination, in such 
cases we instead measured the surface brightness in a 
blank region of sky a few arcseconds away. 

To obtain deprojected galactocentrlc distances of the 
H II regions and B stars, we assumed a constant position 
angle of 22° and a constant inclination angle of 50° (see 
ii2.3l) . The ellipses were centered on the galaxy nucleus at 
0l''33™50!90 -F30°39'36'.'2, which is consistent with the 
hsted position o f M33 in the 2 MASS Extended Source 
Catalog (see also lJarrett et a l."003') . The surface bright- 
ness profile that we derive is displayed in Figure [21 along 
with an exponential disk fit. The large-scale background 
subtraction that was performed on the image makes ac- 
curate assessment of surface brightness errors extremely 
difficult. We therefore added a minimum uncertainty of 
0.1 mag arcsec^^ in quadrature to the measured Poisson 
uncertainties (T. Jarrett 2004, personal communication). 
Our fits to the surface brightness profile (see i)2.3|l yielded 
an extrapolated central surface brightness for the disk of 
liK = 17.51 ± 0.06 and a disk scale length of 5.86' ± 0.28' 
(1.36 kpc). This value for the scale length is in agreement 
with the previous nea r-infrared photometric analysis by 
lR,egan fc Vog-ell l|1994|) , but the central surface brightness 
we measure is ~ 0.15 mag brighter, even after accounting 
for the different inclination angles. Because M33 is very 
extended compared to the largest near-infrared CCD mo- 
saics currently available, determining an accurate pho- 
tometric zeropoint for the surface photometry is quite 
difficult, which is likely the cause of the different central 
surfac e brightnesses derived by us and Rcean & Voeel' 
l|1994D . The imaging by Regan & Vogel (1994) did not 
reach the edge of the galaxy, suggesting that they may 
have derived an incorrect sky level, but the sky subtrac- 
tion is sufficiently problematic that either or both works 
may have a systematic zeropoint offset of ~ 0.15 mag. 

2.3. Rotation Curve 

In order to calculate the dynamical mass of M33, we 
need a rotation curve. The best existing r o tation curves 
of M33 are those constructed by iCorbellil lj2003(l in CO 
and [^prbelli & Salucci (2Q^ in H I, which are both 
consistent with t he hig h-resolution CO observations of 
lEneargiola eTal] l|2003[i . iCorbelli fc Saluccil |2000(l fit a 
tilted-ring model to the H I velocity field that they ob- 
served with the Arecibo telescope. They augmented their 
low resolution (beam FWHM of 3.'5 = 790 pc) data with 
in terferom etric observations near the center of the galaxy 
bv lNewton (198Q ) with a beam size of 47 x 93". The best- 
fitting model constructed by both sets of authors has a 
position angle (PA) of ~ 22° and an inclination angle 
of ^ 50° over the portion of the galaxy covered by the 
2MASS images (both angles change slowly with radius, 
but for simplicity we ignore these variations). Although 
the near-infrared images clearly show that the central 20' 
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Fig. 2.— X-band surface brightness profile of M33 from 2MASS 
data. The error bars represent formal statistical uncertainties from 
the isophotal fits, which clearly overestimate the true uncertainties. 
The solid red line shows the exponential disk fit to the surface 
brightness profile. 



of M33 have a larger inclination angle (58°) and a smaller 
PA (14°), we use the kinematic values of these parame- 
ters for our surface photometry to maintain consistency 
b etwee n the photometry and the mass model. 

iCorb cUi ( 200^ used a new CO map of M33 from the 
14 m FCRAO telescope to confirm the accuracy of the 
H I data and further improve the angular resolution in 
the inner part of the galaxy. We use the CO data out 
to a radius of 6 kpc and the H I data at larger radii. 
The combined CO and H I rotation curve has a resolu- 
tion element that ranges from 45" to 3.'5 and features 
typical velocity uncertainties of ~ 2 km s~^. The over- 
all mass distribution of M33 is very well constrained by 
these data. 

We construct a mass model of the galaxy assuming 
that it consists of a spherical dark matter halo and thin 
disks of gas and stars. The stellar surface density profile 
is derived from the 2MASS Ks image with the IRAF^ 
task ELLIPSE in the STSDAS package, and the atomic 
and molecular sur f ace d ensity profiles are ta ken from 
ICorbelli & Sahlccil ll2000t) and ICorbellil ll200l . respec- 
tively. The rotation curves due to each of the bary- 
onic components are determined directly from the sur- 
face density profiles by numeri cal integr ation using the 
NEMO software package fT\;u ben Ill995|) . The observed 
CO and H I surface density profiles extend out to a radius 
of 80' (the edge of the rotation curve), but the stellar sur- 
face densities can only be measured out to 20'. To extend 
the rotation curve from the stellar disk out to larger radii, 
we extrapolate the stellar surface densities using the ex- 
ponential disk fit from t|2.2l We assume a maximal stel- 
lar disk, which for these data occurs at a stellar mass-to- 
light ratio of M^ILk — 0.51Mq/Lqk- It is worth noting 
that this value is in agreement with the stellar M/L pre- 
dicted from the galaxy color of M^/Lk = 0.55Mq/Lqk 
ijBell fc de Jond |2001j) . The rotational velocities due 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 
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Fig. 3. — Rotation curve of M33. The open circles are the CO 
and H I observations from Corbelli (2003) and CorbeUi & Saluccil 
120001) . The thick gray curve represents the rotation curve of 
the stellar disk (for a stellar mass-to-light ratio of Mt/Lj^ = 
CSIMq/Lq/^ ), the dashed black curve shows the rotation curve 
of the gas disk, and the black rectangles are the rotation curve of 
the dark matter halo after removing the contributions of the stars 
and gas. The thin black solid curve is a power law fit to the dark 
halo rotation curve. 



to the dark halo are then calculated by subtracting in 
quadrature the rotation curves of the stellar and gaseous 
disks from the observed rotation curve (see FigureEl). To 
lessen the influence of small-scale bumps and wiggles and 
ensure that the dark matter rotation curve is monoton- 
ically increasing, we then replace the dark halo rotation 
velocities with a power law fit. We find that Vrot oc r°-^^ 
provides an excellent fit to the data. Since the dark halo 
is assumed to be spherical, the mass profile of the dark 
matter is easily determined from its rotation curve. 

To derive local values of the dynamical M/L, we divide 
each mass component into annuli and add up the total 
mass contained in each annulus. We divide up the light 
profile in the same way, and then take the ratio of the 
mass to the light in each annulus. This M/L profile can 
then be interpolated into a smooth function that gives 
the M/L at any radius. Note that unlike the surface 
brightness measurements, these mass-to-light ratios are 
not truly local measurements in the sense of being deter- 
mined at a single position in the galaxy. Because of the 
nature of a rotation curve derived from a two-dimensional 
velocity field, the M/L values are necessarily azimuthally 
averaged quantities. 

The uncertainties on the mass-to-light ratios that we 
derive deserve some further comment here (see also Q. 
Although the uncertainties on the rotation curve and the 
light profile of the galaxy are known, and therefore we 
can straightforwardly estimate the uncertainty on the in- 
tegrated mass-to-light ratio interior to a particular ra- 
dius, the uncertainty on the local mass-to-light ratio at 
any given radius is less well defined since it is effectively 
the ratio of differences between large numbers. Because 
the other galaxy properties that we are comparing with 
M/L are measured locally, it is important that we use 
local mass-to-light ratios instead of integrated ones, de- 
spite the potentially large uncertainties involved. 
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3. RELATIONSHIPS BETWEEN M/L, METALLICITY, 
AND SURFACE BRIGHTNESS 

In Figures 01 131 and we show that there are strong 
correlations between all three variables under consid- 
eration. We measure correlation coefficients of —0.74 
for the log M/L-metallicity correlation, 0.96 for the log 
M/L-surface brightness correlation, and —0.80 for the 
metallicity-surface brightness correlation. In each case 
the trends go in the expected sense; metallicity and sur- 
face brightness both decline with increasing M/L, and 
surface brightness increases with metallicity. 
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Fig. 4. — M/L-metallicity relationship for H II regions (filled 
squares) and supergiant stars (open squares) in M33. The solid line 
represents the forward fit to the data, the dashed red line represents 
the fit to the inverse relationship, incorporating the metallicity 
uncertainties, and the blue dotted line is the linear bisector of the 
two (Equation Q 
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Fig. 5. — M/L-surface brightness relationship for H II regions 
(filled squares) and supergiant stars (open squares) in M33. The 
solid line represents the fit to the data using M/L as the indepen- 
dent variable, the dashed red line represents the fit to the inverse re- 
lationship, incorporating the surface brightness uncertainties, and 
the blue dotted line is the linear bisector of the two (Equation |2J 



Before we proceed to fitting the data, a brief discussion 
of some of the subtleties involved is warranted. We are 
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Fig. 6. — Surface brightness-metallicity relationship for H II 
regions (filled squares) and supergiant stars (open squares) in M33. 
The solid line represents the fit to the data given in Equation ISl 



in the somewhat unusual situation of investigating rela- 
tionships between three parameters, two of which have 
reasonably well-known and well-understood uncertainties 
(metallicity and surface brightness), and one of which 
does not (mass-to-light ratio). We do not know at this 
point whether one of these parameters physically causes 
the correlations with the others, so we should not assign 
the roles of independent and dependent variables in the 
fits. The most appropriate way of applying the stan- 
dard technique of linear least squa res fittin g is there- 
fore not cl ear. As discussed by Isobe et al.l y.990) and 
lAkritas 'fc Bershad^l ljl996() . there are a number of fitting 
methods that can be used for linear regression problems. 
Because the relationships that we are investigating may 
have intrinsic scatter that is as large as or larger than 
the measurement uncertainties, we use the bivariate cor- 
related e rrors and intrinsic scatter (B CES) method pro- 
posed by lAkritas fc Bershadvl ljl996f) as a generahzation 
of ordinary least squares fitting. In order to avoid bias- 
ing the results by designating one variable as independent 
and one as dependent, for each pair of variables we run 
the BCES regression of F on AT as well as the BCES re- 
gression of X on Y, and then take the linear bisector of 
the two as the best estimate of the true relationship. 

3.1. Two- Parameter Correlations 

We begin our analysis by considering the two- 
parameter correlations between M/L and metallicity, 
M/L and surface brightness, and metallicity and sur- 
face brightness. Following PB02, wc fit the mass-to-light 
ratios at the position of each metallicity measurement 
as a function of the metallicities. Using a BCES fit 
with M/L as the dependent variable, we find the follow- 
ing relationship: log{M/LK) = 10.97 ± 1.73 - (1.27 ± 
0.20)[12 -I- log{0 / H)]. Because the metallicities have 
known uncertainties, while the mass-to-light ratios do 
not, it might make more sense to carry out the inverse fit 
with metallicity as the dependent variable, which yields 
12 + log{0/H) = 8.59 ± 0.04 - (0.46 ± Qm)log{M/ Lk)- 
Converting this equation back into the form of the first 
equation, we find logiM/Lx) = 18.69 ± 3.17 - (2.18 ± 
0.37) [12 -t- log{0/ H)], which differs from the forward fit 



at the ~ 3 fj level in both coefficients. The reduced 
for this fit is 1.35. The forward and inverse fit re- 
sults give a sense of the range of fits that are compatible 
with the data. Although there is a clear correlation be- 
tween metallicity and mass-to-light ratio, the large scat- 
ter makes it difficult to define a precise relationship. Be- 
cause the inverse fit incorporates additional information 
(the metallicity uncertainties), it is probably a more re- 
liable estimate of the true correlation. The linear bisec- 
tor (the average of the ffist and third equations) of the 
metallicity-M/L relationship is 

log{M/LK) = 14.83±1.80-(1.72±0.21)[12-|-Zog(O/ff)]. 

The forward and inverse fits, along with the linear bisec- 
tor, are plotted in Figure 0] 

If we relax our assumption that the stellar disk of M33 
is maximal and construct mass models with lower stel- 
lar mass-to-light ratios, we find only modest changes in 
the derived relationship. Even for the limiting case of a 
minimal disk (in which the stars do not contribute any 
mass to the galaxy), the fit parameters change by at 
most 1.95 (T, in the direction of increasing zero point 
and decreasing (more negative) slope. 

We then carry out the same process with the surface 
brightnesses. For the unweighted fit with M/L as the de- 
pendent variable, we find log{M / Lk) ~ —4.96 ± 0.34-1- 
(0.26±0.02)/iif . The inverse fit (including surface bright- 
ness uncertainties), gives = 19.09 ± 0.10 -|- (3.55 ± 
0.23)log{M / Lk) and when we rearrange this equation 
to write it in terms of M/L, we get log{M / Lk) = 
-5.38 ± 0.34 -I- (0.28 ± 0.02)/iK, this time in reasonably 
good agreement with the forward relationship. The for- 
mal reduced value for this fit is quite high (9.3) be- 
cause of the small errors on the surface photometry near 
the center of the galaxy, indicating that the observed 
scatter is intrinsic to the relationship. Nevertheless, it is 
clear that the fitted lines are a good description of the 
data (see Figure EJ. The linear bisector of the surface 
brightness-M/L correlation is 

log{M/LK) = -5.17 ± 0.24 + (0.27 ± 0.01)/!^. (2) 

For the M/L-surface brightness relationship, making 
the stellar disk less massive does have a statistically sig- 
nificant effect on the derived parameters. The linear bi- 
sector correlation for a submaximal {Vdisk/Vtot = 0.6) 
disk with A'U/Lk = 0.27Mq/Lqk is log{M/LK) = 
—5.76 -I- 0.30^if , and for a minimal stellar disk the cor- 
relation is log{M/LK) = -6.50 + QM^k- 

Finally, for the metallicity-surface brightness relation, 
we obtain a best fit (taking uncertainties in both param- 
eters into account) of 

12 + log{0/H) = 10.99 ± 0.40 - (0.13 ± 0.02)/!^, (3) 

with a reduced value of 1.02 (see FigureEl). Note that 
in this case, since the observed scatter is consistent with 
being entirely attributable to the observational uncer- 
tainties, we do not need to use the BCES technique. We 
instead employ the method of Press ct al. ( 1992) for lin- 
ear fitting with errors in both coordinates. The inverse 
fit gives identical results for this relationship (which is 
not surprising, since the variables are treated completely 



6 



Simon et al. 




-0.5 0.0 0.5 1.0 

-6.62 + 0.15[12 + log(0/H)] + 0.28/^^ 



Fig. 7.— The local fundamental line within M33. The filled 
squares represent H II regions and the open squares represent su- 
pergiants. The solid line shows the fit to the data given in Equation 



symmetrically). This result is in excell ent agr e ement 
with the /-band relationship derived by iRvded ijlOQIif l 
from measurements of 97 H II regions in six galaxies, 

12 + log{0/H) = 10.91 ± 0.90 - (0.095 ± 0.034)/i/. (4) 

3.2. The Fundamental Line 

Having established the various two-parameter correla- 
tions between the data points, we now consider the best- 
fitting plane for all three parameters simultaneously. Us- 
ing an unweighted ordinary least squares fit, we find that 
the fundamental line can be expressed as 

log{M/LK)^-&.Q2±lA2 

-H(0.15±0.16)[12 + ;o.g(O/i/)] (5) 

-|-(0.28±0.03)/iA'. 

The local fundamental line for M33 is plotted in Figure 
d The effect of increasing the assumed distance for M33 
from 800 kpc to 900 kpc (see ^2.1|l is to slightly increase 
the constant offset in Equation |5l to —6.47. The metal- 
licity coefficient is unchanged, and the surface brightness 
coefficient decreases to 0.27; all three terms change by 
much less than 1 a. Adjusting the mass model to allow 
for a lower mass stellar disk also does not significantly 
change the fundamental line. Even for a minimal stel- 
lar disk the constant term and the metallicity coefficient 
change by less than their uncertainties, and the surface 
brightness coefficient increases by 1.7 a. The derived 
fundamental line therefore appears to be robust against 
plausible systematic uncertainties in our assumptions. In 
FigurelHwe display edge-on and end-on views of the fun- 
damental line in three dimensions, proving that the rela- 
tionship is indeed a line and not a plane. 

In contrast with the dwarf galaxy data, we find a very 
weak metallicity dependence for the fundamental line in 
M33. The metallicity dependence is also in the oppo- 
site sense of what we found when considering the two- 
parameter M/L- metallicity relationship. This may be 
partially caused by the small range of metallicities in- 
cluded in our sample (~ 1 order of magnitude) and the 
relatively large uncertainties on those metallicities, but 



could also indicate a real difference between M33 and the 
Local Group dwarfs. In support of the former interpre- 
tation, it is worth noting that the M/L-metallicity cor- 
relation is the weakest of the three relationships studied 
in t)3.1| Measurements of a large sample of metallicities 
including the innermost and outermost regions of M33 
may be able to determine whether the metallicity depen- 
dence of the disk galaxy and dwarf galaxy fundamental 
lines actually differ. 

3.2.1. The Curvature of the Fundamental Line 

Despite the name that we have chosen for the funda- 
mental line, inspection of Figure suggests that rather 
than being truly linear, the fundamental line is actually 
slightly curved. At both the low M/L and high M/L 
ends of the line the data are curving upwards, towards 
larger mass-to-light ratios than the linear relation would 
predict. These deviations from linearity can be naturally 
understood by considering what the local values of M/L 
are physically describing. Near the center of the galaxy 
(at low M/L), the baryons contribute substantially to 
the total mass, increasing the mass-to-light ratio. If we 
used only the dark matter halo to calculate masses, then 
this end of the line would be more nearly straight (as 
we show explicitly in t|4.2|l . The other end of the funda- 
mental line also curves upward because there is almost 
no light being added at large radii, making the mass-to- 
light ratio increase more rapidly. Consistent with these 
explanations, we note that similar curvature is evident 
in the M/L-metallicity and M/L-surface brightness cor- 
relations, but not in the metallicity-surface brightness re- 
lationship, indicating that the origin of the curvature is 
in the mass-to-light ratio. Because of the limited extent 
of the available data, it is more useful to use the linear 
description of the fundamental line that we have derived 
rather than a quadratic one. However, future studies 
that employ larger data sets may be able to examine the 
curvature of the relationship in more detail. 

4. DISCUSSION 

In the previous section, we showed that not only are 
mass-to-light ratio, metallicity, and surface brightness all 
correlated with each other in M33, but also that a local 
fundamental line, which holds at every measured point 
within the disk of the galaxy, exists as well. We now ex- 
amine this relationship in comparison to the dwarf galaxy 
fundamental line and consider its implications. 

4.1. Comparison with the Fundamental Line oj Dwarf 

Calaxies 

Since the motivation for this work was the discovery 
of the fundamental line of dwarf galaxies, we would like 
to compare the relationships that we have found within 
M33 with the analogous relationships among dwarfs. Us- 
ing linear least squares fits to the forward relationships, 
PB02 derived the following three correlations for the in- 
tegrated properties of Local Group dwarfs: 

log{M/Lv) = -QA-[Fe/Hl (6) 

log{M/Lv) = -5.35 + Q.21^iv, (7) 

and 

log{M/Lv) = -2.86 - 0.56[Fe/H] + Qli^iy. (8) 
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Fig. 8. — Edge-on (left) and end-on (right) views of the fundamental line in the three-dimensional space (metallicity, surface brightness, 
log(M/L)). The narrowness of the fundamental line when observed end-on demonstrates that the points are indeed distributed in a line 
rather than a plane. 



If we fit the same data using the BCES method, as we 
did for M33 in ^'6.1\ we find that the coefficients differ at 
the ^ 1 — 2 a level: 

log{M/Lv) = -0.68 ± 0.13 - (1.13 ± Om)[Fe/H] (9) 

and 

log{M/Lv) = -6.23 ± 0.53 + (0.31 ± 0.02)^y (10) 

In order to compare these fits properly with our results, 
we must derive conversions between the [Fe/H] and [12 + 
log(0/H)] abundance scales and V- and iC-band mass- 
to-light ratios and surface brightnesses. A subset of the 
Local Group dwarf galaxies discussed bv iMateol l)1998f ) 
have both stellar [Fe/H] measurements and H II region 
oxygen abundances, from which we obtain 

12 -I- log{0/H) ^ 9.5 ± 0.5 + (1.1 ± 0.4)[i^e/i?], (11) 

or alternatively, 

[Fe/H] = -8.5±2.6+(0.9±0.3)[12-HZo5(O/i7)]. (12) 

The abundance range spanned by the Local Group dwarf 
galaxies encompasses the abundance range of our data in 
M33 (see Figure EJ. 

if-band observations do not exist for many of the Lo- 
cal Group dwarfs, so the V- to if-band conversions must 
be derived specifically from the M3 3 data. M33 has inte- 
grate d magnitudes of my ~ 5.37 ijde Vaucouleurs et al.1 
1199111 and tuk = 2.73 (both measurements are cor- 
re cted for Galactic and interna l extinction according 
to 'Schlcge l. Finkbeiner. fc Davi^ [1998] and Saka i et all 
[2000], respectively), giving a.V — Kg color of 2.64. We 
can therefore estimate that 



\ogM/Lv = \ogM/LK + 1.06 (13) 

and 

f^v = /i/f + 2.64. (14) 

Applying the conversion given in Equation 1121 to the 
dwarf galaxy metallicities, we find that the dwarf galaxy 
M/L- metallicity relation (Equation |5J can be written 

log{M/Lv) = 8.93±3.04-(1.02±0.35)[12 + Zo5(O/i/)]. 

Converting the M33 relationships to V^-band using Equa- 
tions El and El gives M/L- metallicity and M/L-surface 
brightness relations of 

log{M/Lv) = 15.89±1.80-(1.72±0.21)[12 + ;o5(O/i7)] 

(16) 

and 

log{M/Lv) ^ -4.82 ± 0.24 + (0.27 ± 0.01)nv (17) 

Comparing Equations El f^^d El we see that both the 
slope and the zero point of the M/L-metallicity rela- 
tionship are only marginally consistent (at the ~ 2 cr 
level) between the Local Group dwarf galaxies and M33. 
The agreement between the global and local M/L-surface 
brightness relationships (Equations 1 1 01 and 1 1 7|l is slightly 
worse, with a slope difference of 2 cr and a zero point 
offset of 2.7 (7. However, as shown in Figures [51 and El 
for all three correlations the trends in the M33 data do 
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Fig. 9. — Comparison of the M/L-metallicity (left panel) and M/L-surface brightness (right panel) relationships within M33 and among 
the Local Group dwarfs. The M33 data are plotted as filled squares, and the dwarf galaxy points are plotted as open circles. The solid 
lines show the linear bisector fits to the M/L-metallicity and M/L-surface brightness relationships from Figures Bl and O The dashed 
(submaximal disk with M^/Lx = 0.27 Mq j Lqk) a^nd dotted (minimal disk with M^/Lx = Mq / Lqh) lines illustrate the effect of 
using different mass models for M33. Lowering the assumed stellar mass-to-light ratio does not improve the agreement between the M33 
and dwarf galaxy data points. 
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Fig. 10. — Comparison of the surface brightness-metallicity re- 
lationships within M33 and among the Local Group dwarfs. The 
M33 data are plotted as filled squares, and the dwarf galaxy points 
are plotted as open circles. 



Fig. 11. — Comparison of the fundamental line within M33 and 
the fundamental line of the Local Group dwarfs in the same units. 
The M33 data are plotted as filled squares, and the dwarf galaxy 
points are plotted as open circles. 



appear roughly consistent with the trends among the Lo- 
cal Group dwarfs except for a zero-point offset. The local 
M/L in M33 is roughly a factor of 3 higher than the global 
M/L of a dwarf galaxy at the same surface brightness, 
and a factor of 5 — 10 higher at the same mctallicity. It 
is worth noting that using different mass models for M33 
that place less mass in the stellar disk does not improve 
the agreement between the local and global relationships. 

Converting the M33 and dwarf galaxy fundamental 
lines to the common system of oxygen abundances and 
V-band magnitudes yields 



log{M/Lv) ^ 1.90-0. 50[12 + log(O/H)] + 0.13Aiy (18) 
for the dwarf galaxies, and 



log[M/Lv) = -6.30±1.424-(0.15±0.16)[12+log(O/i/)] + (0.28±0.03 

(19) 

for M33. These relationships are clearly in rather poor 
agreement with each other. Nevertheless, a visual com- 
parison again demonstrates that the variation of the 
mass-to-light ratio with mctallicity and surface bright- 
ness is similar in M33 as it is for the dwarf galaxies (see 
Figure [TT|) . 

The similarities in the slopes of all of these relation- 
ships indicate that there could be a connection between 
the local and global fundamental lines. However, M33 
does not lie on the fundamental line of dwarfs, so the 
nature of this connection is not clear. Because struc- 
ture forms hierarchically in the currently-favored Cold 
Dark Matter (CDM) cosmology, disk galaxies such as 
M33 were presumably constructed from the merging of 



Is There a Fundamental Line for Disk Galaxies? 



9 



many smaller galaxies. We suggest that simple simula- 
tions could be used to test whether a galaxy that is built 
up of smaller objects that obeyed the fundamental line 
of dwarf galaxies follows a local fundamental line like the 
one we observe in M33 at the present day. If not, what 
additional information is necessary to create the funda- 
mental line? 

4.2. The Origin of the Fundamental Line Correlations 
4.2.1. The Metallicity- Surface Brightness Relation 
The correlation between local surface brightness and 



ous work (e.g ., Edmunds & Paeel 1984; 



ship is essentially a generic result of ordinary chemical 
evolution models. For example, if we treat M33 as a 
closed box, then the heavy element fraction Z as a func- 
tion of the gas fraction /i can be written as 



L from 


previ- 


Rvder 


1995; 



(20) 



where p is the yield of h eavy elements from eac h gener- 
ation of star formation (|Searle fc Sargeni]ll972fl . If we 
write out the expressions for Z and /x, then we have 



Mz 



-p In 



M. 



gas 



Ma. 



AL 



(21) 



^gas \-^*-^gas 

where Mz is the total mass of heavy elements and M» 
is the stellar mass. Now, we divide the mass terms on 
the right-hand side of the equation by a unit area to turn 
them into surface densities: 



In this case also, for a yield oip = 4.5 x 10^"', we correctly 
recover the surface brightness-metallicity relationship in 
M33. It appears likely, then, that most plausible chemi- 
cal enrichment scenarios will produce a strong correlation 
between metallicity and surface brightness. We therefore 
expect that all disk galaxies should display such a corre- 
lation. Galaxies that lie on the same corr elation, as M33 
and the sample studied bv iRvdeil 1)19951) do, may have 
similar values for the heavy element yield. 

4.2.2. The M/L-Surface Brightness Relation 

Relationships between surface brightness and mass-to- 
light ratio have been recognized before on a global level 
(e.g ., Z waan et al. 1995; Zavala et al. 2003, although 
see iGrahanJ l2002fl , but the corresponding relationship 
withi n a single gal axy has received only minimal atten- 
tion l)Petroulll98H) . In order to understand the origin 
of the local M/L-surface brightness correlation that we 
found in t l3.1l we consider a galaxy wit h an exponential 
stellar disk and a iNavarro. Frenk. fc White! <jl996l here- 
after NFW) dark matter halo: 



I{r)=Ioe-''/'''' 
p[r) = - 



(27) 
(28) 



Integrating the surface brightness profile over the circular 
disk area and the density profile over a spherical volume 
give the integrated luminosity and mass profiles of the 
galaxy, respectively: 



Mz 

Mgas 



= —pin 



-'gas 



-'gas 



(22) 



After substituting the oxygen mass for the total heavy el - 
ement mass, using Mq = 0A5Mz ijGarnett et al.llT99^ . 
and noting that Mgas = ISSMh, we can rearrange Equa- 
tion|211to obtain the following for the metallicity-surface 
brightness relationship: 



^gas 



11.78. (23) 



(24) 



12 + log (O/i/) = log pin 1^1 
Solving this equation for the yield p gives us 

]^Ql2-|-logO/i/-11.78 

^ ln(l -I- S^/Egas) 

and, using the data from ^ we find that for p — 3.5 x 
10^"* this model reproduces the observed metallicity- 
surface brightness correlation. 

Alternatively, we can abandon the closed-box assump- 
tion and all ow the infall of metal-poor gas onto M33. 
In this case, iClavtonl 1)19871) has shown that the heavy 
element fraction is 



^(m) = |hlnM + ln(-ln/i+l)]. 



(25) 



Proceeding through the same steps as above, we find that 



2^gl2+logO/H-ll, 



48 



In 



In 



1 - In 



(26) 



Lir) = 27rIorl 



M{r)=^TT5cPcrl 



In 1 



1 



(29) 



(30) 



Dividing the galaxy into concentric rings of width Ar, 
calculating the incremental mass (AAf = M{r -\- Ar) — 
M(r)) and light (AL = L{r -\- Ar) — L{r)) in each ring 
using Equations 1291 and 28, and then taking the ratio of 
the two yields 



AAf 



ln(l 



-ln(l + ^) + ^- 



r+Ar 
r+Ar+r,, 



/or^e-'-/'-^ 



1 



J- - g-Ar/rrf 
rd 



1 



r+Ar 
rd 

(31) 

for the local mass-to-light ratio. While this expression 
may look formidable, in fact its logarithm scales linearly 
with surface brightness (and radius) as long as r^ 3> r^; 
(which is always true for real galaxies). Using reason- 
able value s for the disk o f M33 from this work and the 
halo from iGorbellil l)2003D (NFW concentration c = 10 
and scale radius r^ = 15 kpc), we plot the M/L-surface 
brightness relation according to Equation 1311 in Figure 
1121 There is reasonable agreement between the modeled 
and observed relationships. If the model had included 
the baryonic contribution to the total mass, which is sig- 
nificant at small radii (high surface brightnesses), then 
the solid line would curve upward at ^ 19, following 
the data. Using other density profiles for the dark mat- 
ter halo, such as a power law (as we assumed in i)2.3|l 
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or a pseudo-isothermal profile, does not change the es- 
sentially linear correlation between logM/L and surface 
brightness, although in both of those cases the relation- 
ship does bend slightly upwards at high surface bright- 
nesses rather than remaining exactly linear. Thus, the 
M/L-surface brightness correlation appears to be a uni- 
versal property of galaxy disks that simply reflects the 
exponential nature of the disk and the increasing dark 
matter fraction as a function of radius. 
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Fig. 12. — Model M/L-surface brightness relationship for a 
galaxy with an NFW dark matter halo an d an exponential disk. 
The solid line shows the log of Eguation lijll plotted against surface 
brightness, using parameters for the dark matter halo and stellar 
disk that are appropriate for M33. The solid points show the ob- 
served data for M33, and the dashed line gives the M/L-surface 
brightness relationship that we derived from the data in Equation 

H 

4.2.3. The M/L-Metallicity Relation 

Although the metallicity-surface brightness and M/L- 
surface brightness relationships follow generically from 
the basic processes that occur in galaxies, what is not 
clear is whether the M/L-metallicity relationship like- 
wise is fundamental to the nature of disk galaxies, or if 
it is merely a residual correlation resulting from the other 
two relationships. Since the M/L-metallicity relationship 
has the largest scatter of the three, one might be justi- 
fied in suspecting the latter. However, as we discussed 
in ^ PB02 proposed a model that can reproduce the 
global M/L-metallicity relationship in dwar fs, and simu- 
lation s confirm that this result is plausible ijTassis et al.l 
l2003|l . Therefore, it may be the case that there is a real 
local M/L-metallicity relationship in disks, but the rela- 
tionship is weakened by the relatively large uncertainties 
on existing metallicity measurements. Future metallicity 
studies should be able to shed light on this possibility. 

4.3. Implications 

Previous work suggested that feedback might be re- 
sponsible for the fu ndamental line of dwarf galaxies. 
IDekel fc Wool lj2003(l studied the scaling relations among 
stellar mass, surface brightness, characteristic velocity, 
and metallicity for both the Local Group dwarfs and 
low-luminosity galaxies in the Sloan Digital Sky Survey 
database. They argued that a simple recipe for super- 
nova feedback could reproduce the dwarf galaxy funda- 



mental line. iTassis et al.l l|2003j) used cosmological hy- 
drodynamic simulations to show that a similar feedback 
model produces a M/L-metallicity relationship for dwarf 
galaxies consistent with that found by PB02. It there- 
fore seems plausible that the same type of model might 
be able to explain the local fundamental line within M33 
(no te that M33 do es lie within the mass range for which 
the IDekel fc Wool scaling relations apply, with a stellar 
mass of 4 X 10^ Mq, compared to the mass of 3 x 10"^° Mq 
where the relationships change form). However, our anal- 
ysis in H4.2l indicates that at least two of the three correla- 
tions that comprise the local fundamental line can be de- 
rived without appealing to feedback. Furthermore, other 
numerical work shows that both the local and global fun- 
damental lines can be produced in simulations that do 
not contain feedback f ii4.4l of this paper and A. Kravtsov 
2005, p ersonal communicati on) . 

If the lDekel fc Wod l)2003|) picture is correct, then only 
galaxies below the transition mass of 3 x 10^" Mq in 
stars (~ 6 X 10"'^^ M0 in total mass) should exhibit a lo- 
cal fundamental line. For more massive systems, Dekell 
(2004) suggests that feedback from active galactic nu- 
clei plays an important role in shaping the evolution of 
galaxies. A critical test of the supernova feedback model 
is therefore whether a fundamental line relationship holds 
for massive, early-type disk galaxies and massive ellipti- 
cals. This test could be carried out for nearby galaxies 
with photometry of large samples of red giant branch 
stars (to measure the mean metallicity as a function of 
radius) along with long-slit spectroscopy to obtain the 
velocity dispersion profile (and hence M/L). 

A further implication of the local fundamental line is 
that if low surface brightness (LSB) disk galaxies fol- 
low the fundamental line, then they should have large 
mass-to-light ratios and low metallicities. Observations 
clearly indicate that LSB galaxies are metal-poor (e.g., 
Ide Narav et al.l 12004^ ). and there is also evidence for 
higher mass-to-light ratios than are found in normal 
galaxies llZwaan et alJ 119951 Ide Blok fc McGaugbl 119961 
although see lSpravberrv et alJll995|) 

Finally, the fundamental line of M33 can also be used 
to place constraints on CDM simulations of the formation 
of disk galaxies. Hydrodynamic simulations of galaxy 
disks are beginning to approach the actual properties 
of observed disks in terms of size, mass, and rotation 
velocity, although t he formation of bti l geless disks such 
as M33 is rare (e.g.. |Abadi et al. 2003VRobcr tson et alJ 
12004: .Gove'rnato et alj|2004) . The results of this paper 
show that realistic simulations of disks must also obey 
the fundamental line; if they do not, then the star for- 
mation and/or feedback recipes they include must not be 
correct. Current simulations are already detailed enough 
to investigate mass-to-light ratios, surface brightnesses, 
and metallicities on small physical scales, so in the fol- 
lowing subsection we carry out this test. 

4.4. The Fundamental Line in Simulations 

We analyzed the fundament al line properties of the 
ACDM disk galaxy simulation of lRobertson et ah) 1 20041 , 
updated using the GADGET2 code (lSDringe]"2005V. The 
simulation uses the Soringcl fc Hcrnauist (2003) model 
for the multiphase ISM to explore the impact of ISM 
pressurization by supernova feedback on disk galaxy for- 
mation, and produces a disk galaxy with a mass of 
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2.16 X lO^^ h^^ M(7i and an expone ntial surface brightness 
profile (see lR.obertson et "311120041 for further details). 

We extracted mass-to-light ratios, K-band surface 
brightnesses, and gas metallicities from the simulation as 
follows. We calculated the angular momentum of the disk 
from the gas in order to define a face-on disk plane. Our 
measurements were then made in the plane of the disk. 
The K-band luminosity of the stellar particles was de- 
termined from the Bruzual & Chariot ( 2003) population 
synthesis models. The metallicity and surface brightness 
were measured in azimuthally averaged annuli, so these 
quantities are not truly local values as our observations 
are. However, the simulated galaxy is fairly axisymmet- 
ric (apart from a weak bar feature), so we do not expect 
this difference to affect the results. The mass-to-light ra- 
tios were calculated by directly measuring the dynamical 
mass in spherical shells and comparing with the luminos- 
ity contained in annuli at the same radii, essentially the 
same method that we used for M33. 

We plot the M/L-metallicity, M/L-surface brightness, 
metallicity-surface brightness, and fundamental line re- 
lationships for the simulated galaxy in Figure E| All of 
the trends we found in the M33 data in |Slare matched by 
the simulation, although there are also clear offsets be- 
tween some of the actual relationships and the simulated 
ones, most notably involving the metallicity. Remark- 
ably, though, the fundamental line is reproduced very 
nearly correctly in the simulation. 

The gas metallicities seen in the simulation are too 
large by ^ 0.6 dex, which is not surprising because 
the simulation does not include a model for strong 
supernova-driven gas outflows. Also, simulated low-mass 
galaxies typically form too many stars given their total 
mass and gas content, thereby producing too many met- 
als as well.^ The zero-point offsets in the left panels 
of Figure El are probably a result of these two effects. 
We interpret these results as an encouraging sign that 
the star formation and feedback models included in re- 
cent galaxy formation simulations require only modest 
adjustments in order to produce realistic stellar popula- 
tions and chemical evolution. 

In order to test whether the feedback in the simulation 
is what causes the fundamental line, we also ran a set of 
simple dissipational collapse simulations of disk galaxy 
formation. These simulations are comprised of 40,000 
gas particles distributed in a fixed NFW dark matter 
halo with a virial mass of Afvir — 2.16 x lO^^h^^ Mq 
(the same as the galaxy in the cosmological simula- 
tion from iRobertson et alJ [2004]). and a concentra- 
tion of 15.24 (appropriate for the mass; IB hillock et al.l 
120011) . The initial galaxy has a gas fraction /gas =0.1 
and is given a spin of A = 0.1. The simulations are 
then evolved for 2 Gyr. In each case, the gas cools 
and quickly settles into a d isk (within ^ 100 Myr) and 
begins to form stars (see iSpringel & HernauistJ 120031: 
IRobertson et all l2004[). W e 'model the ISM in these 
simulations using the Springcl & Hcrnquist (2003) mul- 
tiphase model, modified to include an effective equa- 

® The reason that low-mass galaxies form too many stars in 
simulations may be that the star formation prescriptions are ex- 
trapolated from higher-mass systems. A new formulation based on 
pressure-regulated molecular cloud formation and on observations 
of low-pressure systems produces lower star form ation rates and 
thus lower metallicities IBlitz fc R.osolowskvll2(Kl3) . 



tion of state parameter ((?eos) that linearly interpo- 
lates between an isothermal gas (^eos — 0) a-nd an 
ISM strongly pressurize d by sup e rnova feedback {qeos = 
1) as in the i Robertso' n et al.l ll2004l) simulation (see 
Springcl & Hernauist 2003, Robertson et al. 2004, and 
■Spr i ngcl. Di Matteo. fc H ernauist 2Q0 5 for further de- 
tails). We simulate the dissipational collapse with var- 
ious values of qeos between 0.01 and 1.0 to study the 
impact of the effects of supernova feedback on the fun- 
damental line. 

We find that when gEOS is close to (i.e., feedback is 
weak), all of the scaling relations have significant scat- 
ter (see Figure lCTl . In particular, for the M/L-metallicity 
and surface brightness- metallicity relationships, the scat- 
ter is nearly large enough to wipe out the correlations en- 
tirely. There are also noticeable zero-point offsets in each 
of the relationships, but since these simulations are in- 
tended to be simple models rather than realistic galaxies, 
we do not expect them to reproduce the observed scaling 
relations with high fidelity. Despite the scatter and the 
failure to produce accurate metallicities, a fundamental 
line is still apparent in these simulations. As ^eos is 
increased towards 1, the relationships all tighten consid- 
erably, and the simulated points converge on the real 
data. This experiment suggests that feedback is respon- 
sible for the tightness of the fundamental line, but does 
not set the slope and zero point. It is worth noting that 
if we consider the M/L-metallicity, surface brightness- 
metallicity, and fundamental line relationships in terms 
of the mass-weighted gas-|-stellar metallicity instead of 
simply the ISM metallicity, good correlations (with scat- 
ter) are obtained even in the simulations with very weak 
feedback. In fact, with this metallicity prescription, the 
weak feedback correlations are generally closer to the ob- 
served correlations than the strong feedback results are 
(although again, they contain more scatter). Since the 
gas and stellar metallicities do not agree very well in 
these toy model simulations, it is not clear which metal- 
licity definition offers the best comparison to our obser- 
vations. 

Although the results of these simulations appear 
promising in terms of matching the fundamental line, 
the substantial offsets seen in some of the component 
relationships suggest that we should be cautious in 
our interpretation. In a related example, while it is 
straightforward to generate numerical galaxy models 
that match the slope of the TuUy-Fisher relation, sim- 
ulations have not yet been able to simultaneously re- 
produce t he slope, scatter, and zero point of the rela- 
tion fe.g.. 'Navarro fc Steinniet3l2000t lAbadi et al.ll2003t 
iDutton ct al. 2006). By analogy, it is possible that the 
agreement we have found between the data and the sim- 
ulations is fortuitous. 

5. CONCLUSIONS 

We have shown that there are strong correlations be- 
tween dynamical mass-to-light ratio and metallicity, dy- 
namical mass-to-light ratio and surface brightness, and 
metallicity and surface brightness that appear to hold 
locally, at every point within M33. The trends of these 
relationships are that low metallicity, low surface bright- 
ness regions have high M/L, and high metallicity, high 
surface brightness regions have low M/L. By analogy 
with the fundamental line of dwarf galaxies identified by 
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Fig. 13. — Comparison of the scaling relationships in a simulated galaxy from IRobertson et al.l 120041) with the observed relationships 
in M33. (a) The M/L-metallicity relationship, (b) The M/L-surface brightness relationship, (c) The surface brightness-metallicity 
relationship, (d) The fundamental line. The black squares are the M33 data points, and the green triangles are the measurements of the 
simulated galaxy. All of the trends in the data are correctly reproduced in the simulation, although there is a significant zero-point offset 
in the simulation metallicities. 



PB02, we combine the three correlations into the funda- 
mental line of disk galaxies (Figure Cj), using local prop- 
erties instead of global ones. The M33 data follow the 
same trends as the Local Group dwarf galaxies, but with 
zero-point offsets, such that at a given metallicity, the 
mass-to-light ratio in M33 is approximately a factor of 
5 — 10 higher than that of the dwarfs, and at a given sur- 
face brightness M33 has a mass-to-light ratio that is 3 
times higher than the dwarf galaxies. These conclusions 
are largely independent of the particular mass model cho- 
sen for M33. 

We demonstrated using analytical calculations that the 
metallicity-surface brightness relationship in galaxy disks 
is a generic prediction of ordinary chemical evolution 
models. Similarly, the M/L-surface brightness correla- 
tion can be derived straightforwardly under the assump- 
tions of an exponential stellar disk and a dark matter 
halo that follows any of the proposed structures in the 
literature. It is not yet clear whether the M/L metallic- 
ity relationship shares a similar explanation or is merely 
a residual of the other two correlations. These results 
indicate that at least two of the three relationships that 



make up the local fundamental line originate from the 
basic properties of galaxy disks, and should therefore 
be universally applicable to other galaxies. Feedback, 
which was p r opose d by PB02,|Dekcl & Woo ( 2003), and 
iTassis et al.l l|2003D as the physical mechanism responsi- 
ble for the fundamental line of dwarf galaxies, does not 
seem to play such an important role in the fundamental 
line of disk galaxies. 

Finally, we compared the fundamental line relation- 
ships in M33 with their equivalents in a simulated disk 
galaxy from Robertson et al. (2004). We found good 
agreement in the slopes of each relationship, although 
the metallicity of the simulated galaxy is too high. Addi- 
tional toy model simulations showed that feedback con- 
trols the tightness of the fundamental line, supporting 
the conclusion of our analytical models that feedback is 
not necessary to create the local fundamental line. Forc- 
ing simulated disks to match the fundamental line and 
its scatter will improve current numerical recipes for star 
formation and feedback. 



Is There a Fundamental Line for Disk Galaxies? 
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Fig. 14. — Comparison of the scaling relationships in our dissipational collapse toy model simulations with the observed relationships 
in M33. (a) The M/L-metallicity relationship, (b) The M/L-surface brightness relationship, (c) The surface brightness- metallicity 
relationship, (d) The fundamental line. The black squares are the M33 data points, and the red circles and blue triangles are the 
measurements from the simulations in the cases of weak (gEOS = 0.01) and strong (gEOS = 1) feedback, respectively. The stronger 
feedback model clearly tightens all of the correlations and brings some of them closer to the observed relationships. 
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